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We present a new convergent strong coupling expansion 
for two-level atoms in external periodic fields, free of secular 
terms. As a first application, we show that the coherent de- 
struction of tunnelling is a third-order effect. We also present 
an exact treatment of the high-frequency region, and compare 
it with the theory of averaging. The qualitative frequency 
spectrum of the transition probability amplitude contains an 
effective Rabi frequency. 
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The advent of strong laser pulses has stimulated in- 
terest in strong-coupling expansions in quantum optics 
and quantum electrodynamics. Such expansions are 
also of considerable general conceptual interest in several 
branches of physics. However, particularly in the case of 
periodic and quasi-periodic perturbations, the usual se- 
ries, e.g., the Dyson series, are plagued by secular terms, 
leading to a violation of unitarity when the expansion is 
truncated at any order. In addition, small denominators 
appear in the quasi-periodic case (see the discussion in 
the introduction in |l|). These problems have been for- 
mally solved in a nice letter by W. Scherer [^ and in 
the papers which followed f^Jj^. The main shortcoming 
in these works is that convergence was not controlled, 
an admittedly difficult enterprise. By writing an Ansatz 
in exponential form, and "renormalizing" the exponential 
inductively, we were able to eliminate completely the sec- 
ular terms and to prove convergence in the special case 
of a two-level atom subject to a periodic perturbation, 
described by the Hamiltonian |y 



Hi{t) - ea3-/(i)ai. 



(1) 



The corresponding Schrodinger equation is 

idt^it) = Hi{t)^{t), (2) 

adopting h = 1 for simplicity. Above f{t) is of the form 

fit) = ^F„e™'^*, (3) 

nGZ 

with Fn — F^n, since / is real, and at are the Pauli 
matrices satisfying [cri, (72] = 2i(T3 plus cyclic permuta- 
tions. Assuming Fn of order one, the situation where e 
is "small" characterizes the strong coupling domain. 

It is convenient to perform a time-independent unitary 
rotation of 7r/2 around the 2-axis in p), replacing Hi{t) 

by 



H2{t) = ecTi+f{t)<j3 
and the Schrodinger equation by 

idMt) - H2m{t), 



(4) 



(5) 



with $(i) = exp(-i7rcr2/4)*(i). 

The following result was proved in [|l|. Let / be con- 
tinuously differentiable and g be a particular solution of 
the generalized Riccati equation 



9 -^9 
Then the function $ 

$(i) = 



2ifg + ie^ = 0. 
X ^ C'^ given by 



<t>-it) 



umio), 



where 



U{t) 



R{t){l + igoS{t)) -ieR{t)S{t) \ 



(6) 



(7) 



(8) 



-ieR{t)S{t) R{t){l-ig^S{t)) 



with 50 =5(0), 
R{t) ^ 



exp 



'* / (fir) 
'0 



g{T))dT 



(9) 



and S{t) = J„ R{t) ^ dr, is a solution of the Schrodinger 
equation (^) with initial value $(0) = 



'+(0) 
0-(O) 



A 



simple computation |l| shows that the components 0± of 
^(t) satisfy a complex version of Hill's equation 



+ {±ij' + ^ + f 



= 0. 



(10) 



In [Q we attempted to solve (|T^ using the Ansatz 0(i) = 

exp I —i Jp (/(r) -I- g{T))dT J , from which it follows that g 

has to satisfy the generahzed Riccati equation (|6|). A 
similar idea was used by F. Bloch and A. Siegert in [pi. 

For e = a solution of (||) is given by exp ( ~i L f{T)dT j . 
Thus, in the above Ansatz we are searching for solutions 
in terms of an "effective external field" of the form f + g, 
with g vanishing for e = 0. It is thus natural to pose 



m = E^"G'"^W' 



(11) 



where 



and 



G(")(i) = g(i)c„(i) 



q{t) = exp[i I f{T)dT 



(12) 



(13) 



Inserting (|11|)-(|12[) into (|6|) yields a sequence of recursive 
equations for the coefficients c„, whose solutions are 



ci(i) = aiq{t), 



C2{t) - q{t) 

Cn{t) = q{t) 




i / {ahi-rf - <li.T) '^)dT + a2 



T)Cn^p{T)dT +a„ 



(14) 
(15) 

(16) 



for n > 3, where the a„ are arbitrary integration con- 
stants. The main point is that these constants may be 
chosen inductively in order to cancel the secular terms. 
For instance, in order to cancel the secular term in C2 



in (15), the integrand cannot contain a constant term, 
which equals the mean-value term 



^(^^)-t1-^ 



q^{T)dT ^ 0. 



Then it follows from ( |l5| ) that one must require 
Miafq"^ - q-^) = =^ a? 



M(g2) 



M{q^y 



(17) 



(18) 



It was proved in m that one may proceed in this way 
and establish the absence of secular terms of any order. 
Similar results are valid if (ll^ is not satisfied. 

In the quasi-periodic case we were not able to show 
convergence in ([ll|), and, in fact, it is not expected [Q. 
Hence, (11) is to be viewed as a formal power series. In 
the periodic case (^) much stronger results are possible, 
as we now discuss. 



i(") /-<(") 



(2) 



Let Gm , Cm , Qm and Qm denote the Fourier coef- 



ficients of G(")(t), c„(i) (given in (|llD-(y)), q{t) (given 
in ([l^)) and q^{t), respectively, defined as in (||). Due to 
the multiplication by q{t) in (|lj) the Gm are given by 
convolutions 



r<{n) 



/ ^ Qm — lGi 



(") 



(19) 



i=-c 



.(") 



The Cj , the Fourier components of c„(t), have, by (M)- 

(|l6|), explicit expressions in terms of the Qm and Qm , 
for instance, if (|l^) holds and ai is given by ([iq). 



C(i) 


= aiQm, 






(20) 


^(2) 
^m 




^m~n 


QmQ^-l 


, (21) 



Note that, above, Q^^^ = M{q 
the case Miff') = is found 



in 11 



0. The solution for 
. Finally, let 



n 



n{f) = Fo + Go(e), 



Hn — 



[Fn + Gn{e)]{ncu)-\ 
0, 



for n ^ 
for n = 



(22) 



(23) 



and 7/(e) =iJ2 Hm, with G„(e) = J^ ^in^"- Then 

mGZ n=l 

-R(i), in (i), is given by 



i?(t) 



— p-''7/(«) 



-iOt 



TZ{oo, t), 



(24) 



with 7^(Ar, i) = exp f - ^ i/„e'"'^* j. By (|), the 

\ n=-Ar / 

complete wave function is known once (p3) is given; 
(P) and (p3) also show that the wave-function is of the 
Floquet form, with secular frequencies ifi. 

In reference [[7| we have proven the following result: for 
f periodic the e- expansion ^1j\ ) has a nonzero radius of 
convergence. Our estimate for this radius is not optimal 
and we refrain from quoting it here, but we remark that 
the expansion does converge for high frequencies, i.e., 
w 3> e, a condition that we assume in the following. 

We now consider in (0) the special case 



1 



Fn - ^{Sn.l 



Sn, -l), 

corresponding to f{t) = cos(ajt). For this case 
Qm = Jm I - I and QI^J = Jmi — 



OJ 



(25) 



(26) 



By dlSD, ai = 1. 

The transition amplitude A21 from the lowest energy 
atomic state (") of (|l|) to the upper level (i) is 

A2i{t) = (V'l, C/(i)V'2), (27) 

where U{t) is given by (||) and -02 — -j^ [ li)-, V'l = 
4= ( J ) are the corresponding eigenstates of the rotated 



Hamiltonian H2 , given by 



corresponds to the transition probability 

Miit) = Ui, U{t)i;2 



The tunnelling amplitude 
(28) 



where '(/'2 = (°), V'l = (J). The latter represent the lo- 
calized states (eigenstates of the field term, proportional 
to (T3, in (0)), and A21 — Q means absence of tunnelling 
between these states. Indeed, (0) is a semi-classical ap- 
proximation to the spin-Boson system treated in B. In 
the full quantized case, considered in H, V'l and '02 dif- 
fer macroscopically because they are dressed by photon 



clouds, and for e sufficiently small there is always local- 
ization, i.e., no tunnelling. This is not the case here, as 
we shall see. 

Fig. |l] shows the exact result for 1^21 (i)P to fifth order 
in e for about 12 cycles of to. We see clearly the domina- 
tion of the external frequency w, in agreement with the 
theory of averaging. Eq. (0) may be transformed to 



ith ?/; = exp I i L f{T)dT ) -0, and 



with 



/(ij V') — ( -2ism{ujt)/. 



p2z sin{ujt) / uj 





(29) 



^. (30) 



By ( P9| ) and (pOJ), the averaged equation dt'i'o — e/^°-'('0o) 
with /("HV^o) = T /o^ /(i, V^o) and T = 2^/a;, is 



idtil^o = eJo ( - ) cTiV'o 



(31) 



and a well known theorem 



yields 



^{t) - Mt) 



0{e/uj) on the time scale 1/e. Hence A21 is close to 



^pi, exp I -leJo [—jaitj exp I i as ) -02 



lexp 



ieJn { — ] t ] sin 



sin(a;t) 



Since sin 



sin(wi) 



2 2^ ^2/c4 



fc=0 



sin[(2fc + l)wi]. 



we see that in this case the spectrum is dominated by 
the harmonics of the frequency lu of the external field, 
in agreement with Fig. |[ Notice, however, that, while 
averaging is applicable to times up to 0(l/e), the exact 
theory is applicable to all times. Applying the averaging 
theory to A21, we are led to the matrix element 



sin(ti;i) t / n \ , 

exp I z 0-3 Vi, exp (-teJo (x) (Jit) V^2 



- sin(i^t) 

-je" - sin (e Jo (x) i) 



-« Jo ( - ) sin (e Jo (x) t) , 



with X = 2/a;. This result agrees with ( p4|) and Fig. H, 
which shows the exact result for |A2i(t)|^ to fifth order 
in e for t from to 2n/fl. Fig. || shows that ri(e), the 
secular frequency given by (|22|), dominates in this case. 
There, 51(e) ~ 7.6 10~^ for the values of e and w chosen. 
Notice that by (|l|), (|20l), (H) and ^) n{e) ~ eJo(x), 
to first order in e. Thus, the first order contribution 
approaches zero if x approaches one of the zeros of the 
Bessel function Jo. The second order contribution to 51 is 



e^ y^Q^iCl , and is identically zero, as one sees using 



(pi|). The third order contribution to 51 is 



2e3 



E 



ni , n2— — 00 



J2ni+l{x)J2n2+l{x)J-2{ni+n2 + l){X) 

(2ni + l)(2n2 + 1) 



and is non-zero if x coincides with one of the zeros of 
Jo. Hence, when 2/uj approaches one of the zeros of the 
Bessel function Jq the lowest non-vanishing contribution 
to 51 is of third order in e and, hence, rather small. This 
means that for such values of co the tunnelling is very 
heavily, although not exactly, suppressed. 

Hampering and destruction of tunnelling have been 
studied in ||l^,0 for particles, and in ||l^,|l^ for spins. 
The latter use the method of averaging, but we empha- 
size that in the case treated above, a; ^ e is satisfied 
for e sufficiently small, and thus the result is exact, i.e., 
valid for all times. In addition, the features regarding the 
order of the expansion are new. 

At resonance cu — 2e we are not able to prove that the 
expansion converges. It is, nevertheless, a well-defined 
formal expansion, in contrast to strong-coupling approx- 
imations of Keldish type, which are beset with difficulties 
(see, e.g. ||lj] and references given there). Moreover, as 
we shall show, it includes interesting effects of dressing of 
the atoms by the photon field (in the semi-classical ap- 
proximation) which yields the external field Floquet de- 
scription, rigorously justified in |^. Such effects appear 
in the rotating- wave-approximation (RWA) in the form of 
a Rabi frequency (see, e.g. [|5[), but the present model is 
not close to RWA, since the rotating and counter-rotating 
terms in (|l]) are of the same order of magnitude. More- 
over, the RWA is not justified for large coupling, but the 
solution of (P) might have some similarity to the solution 
obtained when the RWA is performed. If so, the effec- 
tive frequency of oscillation of ^21 would not differ much 
from the Rabi frequency (see [nsj) 



51/ 



[(^ 



2e) 



-4] 



1/2 



(32) 



for w ~ 2e (or, in general, for to — 0(e) and e small). 
Indeed, 51 ^^ makes its appearance in (p4| ) in a most inter- 
esting way: by (|l9|)-(|^)-(|2|) and (|^ iJ„ in 7^(cx), t) 
equals, to first order in e. 



TILJ 



yn(^) 



(33) 



with Fn given by (^51). The greatest contributions of 
(p3| ) arises for small n (due to the factor n~^) and when 
the argument of the Bessel function equals its order, i.e., 
n = 2/10, and the corresponding frequency in ( p4|) is 
nuj — ^uj — 2, which compares well with (P4). In Fig. 
g wc show this last effect for IZi^N, t). We considered 



the quantity £{N, t) 



TZ{N. t) 
7?,(oo, t) 



-1 



which measures the 



error made by including in ( p4| ) only the first N terms of 
the sum involving iJ„ in TZ{oo, t). In Fig. Swe considered 
the resonant case with oj — 2e — 2 10^^ and t — O.Ttt/H. 
The qualitative behaviour is the same for other values 
of t. We see from Fig. that mainly only small n and 
n around 2/cj contribute. The effect of adding in the 
second order contribution is negligible for the range of 
values of e considered. 

In conclusion, the new strong-coupling expansion al- 
lows considerable insight into both the high-frequency 
and resonance regimes, and yields an interesting unex- 
pected result for the coherent destruction of tunnelling. 
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FIG. 1. The amplitude |yl2i(t)|^ Here e : 
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FIG. 2. The amplitude |A2i(t)|^ Here e = 0.01 and uj = 2. 
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FIG. 3. The quantity £{N, t) as a function of A^. 



